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Abstract: 

In order for Active Debris Removal to be accomplished, it is critically important to understand the probable 
rotation states of orbiting, spent rocket bodies (RBs). However, rotational dynamics is non-intuitive and 
misconceptions are common. Determinations of rotation and precession rates have been published that are 
physically impossible. In a state of free precession, the total angular momentum of the object is constant, while 
kinetic energy decreases due to internal friction, approaching rotation about the axis of maximum inertia. For 
solid internal friction the timescale is hundreds to thousands of years for quality factors of -100 and assuming 
metallic rigidities, but for turbulent friction in partially-filled liquid fuel tanks we predict that the preferred 
rotational state is approached rapidly, within days to months. However, history has shown that theoretical 
predictions of the timescale have been notoriously inaccurate. In free precession, the 3-1-3 Euler angle rates (p 
(precession rate of long axis about fixed angular momentum with cone angle 0) and if/ (roll rate around long 
axis) have comparable magnitudes until very close to 0=jt/2, so that otherwise the true “rotation period’’ is not 
simply twice the primary light curve period. Furthermore 0 , nonzero due to friction, becomes asymptotically 
smaller as 0=ti/2 is approached, so that 0 can linger within several degrees of flat spin for a relatively long time. 
Such a condition is likely common, and cannot be distinguished from the “wobble” of a cylinder with a skewed 
inertia tensor unless the RB has non-axisymmetric reflectivity characteristics. For an RB of known dimensions, 
a given value of 0 fixes the relative values of <p and \jr . In forced precession, the angular momentum precesses 
about a symmetry axis defined by the relevant torque. However, in LEO, only gravity gradient and magnetic 
induction torques are dominant, and these cannot cause precession periods shorter than a week, or more likely, 
months. Thus forced precession is probably not observable over observation campaigns spanning a few days or 
less. Building on Smith [7] an expression for the magnetic induction torque on an arbitrarily rotating cylinder, 
hitherto not available in the literature, is presented here. Numerical integrations of the equations of motion for a 
cylindrical RB in LEO with arbitrary initial conditions and subject to magnetic and gravity gradient torques as 
well as prescribed internal dissipation are in progress. Acknowledgements: This work was produced under 
NASA contract NNJ05HI05. 


1. Introduction 

The task of Active Debris Removal [1] will be greatly aided by a detailed knowledge of the rotational states of 
rocket bodies (RBs) which are designated for removal from orbit. For this reason, there has been much recent work 
in the space surveillance community on the problem of determining RB rotation states from light curves (e.g. [2]- 
[4]). However, these rotation states are at present inadequately understood. We focus here on the abundant Soviet 
SL-8’s. The primary phenomena influencing the rotational evolution of RB’s in LEO are the torques due to 
magnetically-induced eddy currents and the gravitational gradient, as well as the effects of internal friction within 
the body. In LEO, torque due to solar radiation pressure is probably at least two orders of magnitude or more 
smaller than the smallest of these effects. 

Much groundbreaking work on the rotational evolution of RBs was performed five decades ago, at the start of the 
U.S. space program. For example, analytical predictions were performed [5,6,7] of rotational evolution of 
spherical and cylindrical RBs under the influence of the eddy current torques resulting from motion of the satellite 
through the earth’s magnetic field. While these analytical calculations are elegant, they are limited to rotation states 
in which the angular velocity remains constant over an orbit period, so magnetic field components can be averaged 
over that period. Furthermore, the torques were only calculated in [7] for the states of spin-stabilization about the 
long axis, and for pure end-over-end (flat spin) tumble. Because of these limiting assumptions, the full spin 
evolution of an arbitrarily-rotating cylindrical RB cannot be investigated with these results, and they certainly 
cannot apply to bodies precessing significantly within an orbit period. 

Torques due to atmospheric drag forces are only comparable to eddy current torques at altitudes below about 160km 
[8], so they are ignorable in this work. The dynamics of satellite orientation under the influence of gravity gradient 



torque is of paramount importance in satellite control, so this influence was investigated in detail many decades ago 
as well (cf. [9]). However, its relevance to interpretation of light curves from RBs in LEO bears investigation, since 
its magnitude can be as large as that of magnetic induction, and that spin-orbit resonances may be a factor in light 
curve interpretation. 


2. “Much confusion about precession” , or “Free and forced precession” 


To understand the rotation states of RBs and their corresponding light curves, it is valuable to examine the essential 
similarities and differences between free precession and forced precession, as well as realistic rates that each might 
exhibit. Since RBs 
typically have nearly 
cylindrical symmetry, 
their rotational dynamics 
are well approximated 
by those of a symmetric 
top, and are treated as 
such here. In a state of 
free precession the 
angular momentum L of 
the body remains 
constant in inertial 
space, while the angular 
velocity to and body 
symmetry axis x 3b 

precess about L at a rate n . ... 

F . f B: inertial space 

commonly referred to as 

the free precession rate: Fig.l: Free precession viewed in the body coordinate system (A) and in inertial space (B). 

The view in (B) is directly relevant to telescopic observations, while (A) is not. 
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where I 0 is the maximum rotational inertia and A I is the difference between the maximum and axial inertia 
components : A/ = I 0 - 1 3 . Thus for positive A/, the precession of to and L about x 3b is clockwise as seen in Fig. la. 
Because of the ease with which this result is obtained from the Euler equations, this is commonly the extent of the 
treatment of free precession in introductory dynamics textbooks (cf. [10]), but it is not useful for observations of the 
phenomenon from inertial space as in telescopic observations. However, if we define a 3-1-3 Euler angle sequence 
using the standard angles tp, 0 and \p (where the inertial z-axis is parallel to L), then by equating the components of 

to in the body coordinate system, found in the above-mentioned development, to the the vector ('p + y/ , also 


expressed in the body coordinate system, the following results are readily obtained: 
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where 
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The roll-rate about the symmetry axis, \j/ , is thus seen to be identical to the free precession rate co 0 , so if !// is 
positive in inertial space (counterclockwise axial roll as seen from above x 3b ), then, for the usual case: 0 < 0 < it/2, 
the x 3b axis also precesses in the positive (counterclockwise) sense about L in the inertial coordinate system, at the 
rate (p as seen in Fig. lb. Both (ft and y/ can be referred to as “precession rates”, leading to the need for 


researchers to explain which meaning they are employing when they refer to free precession. For SL-8s, AI/I 3 -4.3, 
and Eqs (2) and (3) indicate that y/ can be a larger component of the angular velocity than <p , unless 0 is very close 
to jt/2. The above equations can be cast in an even more useful form by using the expressions for kinetic energy E 
and the square of L in the body coordinate system: 


E = 
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and 
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We also define the difference between the actual kinetic energy and its minimum possible value as A E: 


A E = E 
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Combining Eqs. (5) - (7) with (2) - (4) yields 
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Figure 2a: dependence of \j/ and COS 6 on A E . Both quantities maintain substantial magnitudes until very close to 
minimum energy. Note in (b) that even with only 1% of maximum A E remaining, 0 is still 16 degrees from a flat spin, 
and the roll-rate if/ is 45% of ('p . 


Eqs (2), (8) and (9) provide considerable insight into the nature of free precession as the energy is reduced by 
internal friction toward its minimum value at A E =0. Fig. 2a shows the square-root dependence of if/ and cos0 on 
the kinetic energy. Consequently, both quantities remain relatively large until A E is very close to zero. Furthermore, 


it is shown below that the rate of internal dissipation is also a decreasing function of 0, so that loss of kinetic energy 
is probably much slower as 0 approaches jt/2, contributing further to the persistence of a relatively large roll -rate. 

A rather remarkable fact revealed by these equations is that remains constant as energy is reduced. Note that 
since cp is the rate of end-over-end rotation when 0 is close to jt/2, in that common situation it is the angular rate that 
most strongly modulates an RB light curve, so it is frequently identified as twice the “rotation rate”. However, 
since if/ is of comparable magnitude to (p until very close to 0=tc/2, it is incorrect to identify <p with the rotation 
rate. In fact, the true rotation rate can be written 


co = p 2 + if / 2 + 2 <py/ cos# = <p \ 1 + 
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Figure 4:(a) maximum energy state corresponds to axial spin stabilization, (b) minimum energy state is pure end-over-end tumble, 
(c) evolution of if/ versus (p . 

Fig. 3a shows the case of maximum rotational energy (AE=AE max , 0=0) which corresponds to spin stabilization about 
the axis of symmetry, and Fig. 3b shows the case of minimum rotational energy (A.E=0, 0=Jt/2) which corresponds to 
pure end-over-end tumble. Note that in the limiting case (a) when 0 =0, if/ and (f) both represent rotation about the 
symmetry axis, and if/ = (AI/I 3 ) (p . Fig. 3c shows the allowed values of if/ and <p . 

In free precession, L must remain along the inertial z axis. Therefore, horizontal components of angular momenta 
about long and short axes must cancel. This is why if/ and cp have similar values until pure tumble is approached. 

It is easy to show that this requirement leads to recovery of Eq. (3). The practical consequence is that claims of 
large light -curve-derived values of “rotation rate”, with 0 not near jt/2 together with relatively small “precession 
rate” (e.g. [4]) are impossible in a state of free precession. In fact the value of if/ is fixed once cp and 0 are 
determined. 

3.1nternal friction 

The time evolution of the rotation state can be understood through Eqs. (8) and (9) if the rate of internal friction is 
known. While it is difficult to ascribe a precise mathematical expression to the dissipation, the amount and 
mathematical form of dissipation due to various phenomena can be roughly estimated. 


Energy loss due to solid internal friction can be estimated as follows: The time-varying stress within a freely 
precessing cylindrical shell of radius R and length H is due to the difference in centripetal acceleration between 
mass elements close to the angular momentum axis and those far from it, experienced across the diameter of the 
RB. This difference scales as 

Aa c ~ ^r2/?cos<9 (11) 

This leads to a difference in centrifugal force felt by one side of the RB relative to the other of roughly half the mass 
m of the RB times A a , and this differential force is applied across a characteristic area A - Hd (cylinder is hollow) 

, where H and d are the length and thickness of the cylindrical shell. Dividing (m/2) A Cl c by A gives the 
characteristic time-varying stress cr , which varies roughly sinusoidally with the free -precession rate: 


cr = cr 0 cos <X) Q t 


( 12 ) 


where 
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and p is the material density of the cylinder. The accompanying strain scales as 
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Where p is the rigidity of the shell and 1/Q is the phase lag between stress and strain. Q is the quality factor for the 
material of which the shell is constructed. The average value of the stress times the strain rate, integrated over one 
precession period, yields an estimate of the volumetric dissipation rate within the RB due to solid friction. 
Multiplying this rate by the characteristic volume of the RB shell, V rb ~2jrRHd gives the estimate of total dissipation 
rate due to solid friction. This is related to the rate of change of 0 by the chain rule: 

Using Eqs. 2 and 3, the result can then be written 
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i e is ~ 10 2 - 10 3 years for typical parameters relevant to SL-8s and using Q-100 and the rigidityof Aluminum. 
Eq.(15) is easily integrated, giving: 


r 
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The result is plotted in Fig.4a. 
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A more significant source of friction is likely generated by the motion of residual fuel and oxidizer in on-board 
propellant tanks. This is an effect that is of great importance in the aerospace industry, yet analytic attempts to 
quantify its magnitude have met with modest success only when combined with laboratory experiments [11]. It is 
helpful, however, to examine its functional form with scaling arguments. We initially consider the liquid to reside 
within the end(s) of the hollow cylindrical shell, filling the cylinder to a fractional height/. Motion is initially 
assumed to be confined to an Ekman boundary layer of thickness 


A = 



(18) 



where v is the kinematic viscosity of the fluid. For dimethyl hydrazine v = 6xl0" 6 m 2 /sec (Wikipedia reference). 

In the above expression the rotation rate is approximated by (jp . The characteristic amplitude of the fluid velocity 
relative to the interior surface is 


V r „«(2R)l(2nlw) (19) 

The viscous stress in the boundary layer is then approximately 

V~rjV r eilq (20) 


and the strain rate is simply the stress divided by the dynamic viscosity, 
the boundary layer 27tRHS gives an estimate of the total dissipation rate. 


Summing the product over the volume of 
Again using the chain rule, we obtain 
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where p f is the density of the fuel or oxidizer and the fluid motion can be shown to be nonturbulent (Reynolds 
number ~40). i v is a few hours for f~0.1 and using the density and viscosity of dimethyl hydrazine. It has indeed 
been found that transition to flat spin can happen very quickly (within 5 minutes) in an initially spin-stabilized RB 
loaded with a large volume of fuel [11]. However, for the case of SL-8s, the fuel is in auxiliary tanks that are a 
fraction of — 1/8 of the cylinder radius. Since x v scales as the inverse cube of the radius, we may estimate that the 
timescale is ~8 3 times longer than that given by equation x, thus increasing the timescale to perhaps a few months. 

Eq.(21) is readily integrated analytically, giving 

0 = cos 1 [exp (-1 lr v )\ (23) 

Figs. 4a, b illustrate the behavior of the coning angle 0 as functions of scaled time for solid and viscous internal 
friction. While these functional forms should not be assumed to be highly accurate, they serve to illustrate how 0 
may linger for an extended time within several degrees of 7t/2. Consequently, from Eq.(3), a significant axial roll 
rate if/ will also persist. 



4. Free precession versus offset inertia matrix 

Fig. 6: Evolution of 0 with solid internal friction (a) and fluid internal friction (b). Note that 0 may linger close to 90 
degrees for an extended period in both cases. 


It has become evident from observations that a significant number of RBs are not in a rotation state of pure end- 
over-end tumble. The results described above suggest that this may be due to values of 0 that remain significantly 
different from ji/ 2 for extended periods of time. Another possibility that has been recently considered [12] is that 
RBs may have principal axes of inertia that are offset from the symmetry axes by a small angle, due to asymmetric 
mass distributions within the body. An interesting question is whether such motion, defined here as “wobble”, is 
distinguishable from free precession. We simulate wobbling bodies that are in a state of minimum energy rotation 
with a prescribed angular velocity direction as shown in Fig. 5. We define the “home” position of an RB such that 
the (xi b ,x 2 b,X 3 b)- body axes are coincident with those of a chosen inertial coordinate system (see Section 8). We then 
reorient the body through the series of four rotations shown in the figure. The first three rotations bring the RB to a 
state where the y b -axis is directed according to the chosen angular velocity. The fourth rotation, represented by 5 in 
the figure, is spin about the chosen axis. A fifth rotation tips the object by the “wobble angle” referred to as a w 







Fig. 7: Series of rotations used to place an RB in a state of "wobble" - i.e., rotating stabley about the 
maximum axis of inertia, where that axis is offset from a symmetry axis. The direction of the angular velocity 
is specified in inertial space by polar angles X m and (p co . 

These rotations are employed in section 9 to simulate light curves from wobbling versus non-wobbling RBs. 
Examination of Eqs (1) — (3) quickly reveals that the only observational difference between the light curve of a freely 
precessing cylinder with coning angle 0, and a “wobbling” cylinder with a wobble angle of tt/2-0. will be if the 
cylinder has a non-axisymmetric reflectivity pattern. Further research will investigate whether the larger cross- 
sectional area when saddle tanks of an SL-8 are visible, leads to enough of a modulation of the signal intensity to be 
readily identified. Initial estimates indicate an approximately -25% increase in reflected intensity due to that effect. 
A useful fact is that if (p and 0 are determined from the light curve data, the signature of the axial roll should have 
the prescribed frequency given by Eq.(3). This signature could of course be sought in FFT data of light curves. 
Another possible signature that would identify a non-zero roll rate would be the insignias on one side of the RB. For 
protection, such insignias are emplaced on the side opposite the launch tower [13] but it is not known how long they 
persist under the harsh conditions of launch and residence in space. 


5. Characteristics of forced precession 


In a state of forced precession, the L 
vector of a rotating body executes a cone 
about an axis in space that is determined 
by the nature of an applied torque. This 
axis forms the natural direction for the 
inertial -axis in the standard 3-1-3 Euler 
angle sequence, so that (ft is along that 
axis. Note that this involves a different 
choice of coordinate systems that of free 
precession, because here the z-axis is 
defined as the axis about which L 
precesses, rather than the L-axis itself. In 
the standard analysis of a symmetric top 
undergoing forced precession, t/7 is 
directed along the symmetry axis of the 
top (as in free precession), it represents 
the largest component of spin, and it is 
approximately equal to the rotation rate. L 
therefore points nearly along the symmetry 
axis (Fig. 6a). However, for forced 

precession of a tumbling RB, the standard Euler angles are not particularly convenient, since precession is about a 
transverse axis so axial roll must be treated by varying 0. 



Figure 8: Standard Euler angle development for forced precession of 
symmetric top (a) is not convenient for the typical case of a tumbling RB. 


The rate at which L precesses about the symmetry axis due to torque t is given by 

d<j> _ r 

dt Lsin# ^4) 

There are only two torques that are relevant for tumbling RBs in LEO: those due to magnetically induced eddy 
currents, and the gravity gradient torque. These are treated in the following sections. 


6. Gravity gradient torque 


The gravity gradient torque exists due to the slight change in the force of gravity across the length of satellite. For 
our case of a symmetric cylinder, the magnitude of this torque can be expressed as: 
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where a g is the angle between the long axis and the geocentric vector, defined here as the unit vector directed from 
the center of mass of the cylinder toward the earth’s center. The direction of the torque is given by the cross product 
of the x 3b body axis and the geocentric vector. With point masses m i and m 2 placed at the ends (see Fig. 8a), the 
above result is replaced by 


A/„=A/ +Mz 2 +m, 


(26) 


where 
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This modification allows investigation of the effects of adding masses to the cylinder. 

For a satellite in circular orbit, it is easily shown that the equation of motion for the angle 2a„ is that of a pendulum 
in the angle 2a g : 


/ 0 «, 
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If the angular velocity is close to the orbit period, Eq. (27) gives simple harmonic motion about either a„ 
with frequency: 

1 2GM A1 

“W V 3 

For an SL-8 at an altitude of 600km, the period of these oscillations is 

P= — «l.3hr 

The period of oscillations is an excellent, intuitively useful measure of the strength of the torque: gravity gradient 
torque only becomes important for SL-8s when rotation period slows to -1.3 hours, when capture into the state of 
synchronous rotation is imminent. Note that this period is very close to the orbital period in LEO, suggesting that 
spin-orbit resonances may be possible for such satellites. Interestingly, the ratio of gravity gradient oscillation 
period to the orbital period depends only on the satellite inertia components: 


= 0 or 7t 
(29) 


f 7 0 V -0 O for SI .-8s (31) 

P 0 m, Ua/J 


It may be that slowly rotating satellites in orbits with nonzero eccentricity exhibit large 
oscillations about either of the stability points (a g = 0,ji) due to the resonant 
interactions with the time-varying gravitational force. Given the estimates of 
components for SL-8s used in this paper, it also appears that resonant 
commensurabilities exist near ratios of orbit period to rotation period of 4/5 and 7/8. 
The practical consequence of these interactions would be that when rotation periods 
are close to the orbit period, chaotic rotation or large librations about these points may 
occur. We should not be surprised to find that the symmetry axes of slowly rotating 
RBs are directed away from the geocentric direction. 

Gravity gradient torque may also produce a slow precession of the rotation axis of a 
more rapidly rotating satellite (rotation period « orbit period) that is nearly in a 
state of flat spin. The situation is shown in Fig. 7. Using the results described 
previously, the precession period is: 



direction 


Fig.9: Precession geometry due to 
gravity gradient torque acting on 
tumbling RB. 
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For an SL-8 with a rotation period of 60 seconds at an altitude of 600km, the precession period is nearly two weeks. 
For typical RBs, it is clear that precession rates due to gravity gradient torque is unlikely to generate precession rates 
that are observable within observations made during a single orbital flyby. 


7. Magnetic eddy current torque 

Electrical eddy currents are induced in rotating, electrically conducting satellites through Faraday’s Law, as the flux 
due to the geomagnetic field passing through the object continually changes. Through the Lorenz force, these 
currents result in a torque on the satellite that is a linear function of the angular velocity components. Components 
of the torque oppose the angular velocity while others are at right angles to the components of omega. Therefore 
this torque causes damping of omega and concommitant precession. Eddy current torque depends on the cube of the 
geocentric distance, so they are most significant in LEO. Numerous papers have been written on the topic of eddy 
current torques acting on cylindrical satellites, but all assume a specific rotation state, and furthermore the torques 
are typically averaged over the orbit period, which limits the results to cases where the angular velocity changes 
slowly relative to the orbit period. Smith [7] computed the torques for two specific cases: axial spin and pure tumble. 
Due to the linearity of the equations, careful examination of the differing coordinate systems and variables defined 
by Smith allows a single expression to be deduced, that is valid for all orientations of the angular velocity relative to 
the RB. In this manner, we find: 
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In this equation the Bi are the components of the magnetic field, measured in the body coordinate system. The 
constants in A’ are: 


ct = electrical conductivity of shell 
t = thickness of cylinder walls 
1 ,R= length and radius of rocket body 

The precession generated by the magnetic torque has a period given by 
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For the same SL-8 parameters described in the previous section, and given the earth’s dipole moment is m~8xl0 22 
Am 2 , P mag is on the order of a few months, and so it is comparable to the gravity gradient torque. Due to this torque, 
it is found in [7] that typical satellite spin components in the equatorial plane of the earth damp in -9 days, and ~90 
days for spin axis parallel to the earth’s dipole axis. 


It is of academic interest to note that x mag goes to zero only when the RB has zero angular velocity, yet the 
gravitational torque vanishes only for a finite (synchronous) rotation state. It can be shown that an RB in a simple 
polar orbit will consequently have a stable state in which a g is slightly nonzero (-.01 degrees), while a negligibly 
small but finite amount of orbital energy is continuously transformed into heat via the eddy currents, causing the 
orbit to experience very slight decay. 

8. A coordinate system framework for simulations of rotation states 

In order to perform numerical integrations of the equations of motion and to relate the rotational dynamics of an RB 
to telescopic observations, a rotational evolution code was created in Matlab. The code includes the orbital 
elements, observation site coordinates, the Greenwich Mean Sidereal Time of the observations, and a set of 
coordinate systems that allow transformation of vectors between the body coordinate system, the coordinate system 
of the earth’s dipole field and of the observation site, and other intermediate coordinate systems. Space does not 
allow a full description here, but by defining a few important coordinate systems, it will be possible for the reader to 
understand our preliminary results. The body coordinate system is defined such that the x3b-axis is directed along 
the symmetry axis of the RB, and the other two axes correspond to the two maximum principal axes (Fig. 8a). The 
inertial “I” coordinate system, (x : , y : , z : ) is defined such that the x : and y t axes are in the equatorial plane, with the 
x r axis directed toward the satellite’s ascending node (Fig. 8b). 

Also shown in Fig. 8b is the gamma coordinate system, where the x and y axes are in the earth’s equatorial plane, 
with the x-axis directed toward the vernal equinox. In order to include the magnetic eddy current torque in the 
equations of motion, the position vector of the RB is transformed to a coordinate system in which the earth’s 
magnetic dipole field B is symmetric about the z-axis (not shown). The components of B are then computed at the 
satellite’s location, and subsequently transformed to the body coordinate system, where they are input into the Euler 
equations. Numerical integrations are in progress. Figs. 9a and 9b illustrate simulated light curves from an open- 
ended cylinder rotating in its minimum energy state (maximum axis rotation) , but here we show only the light 
curves corresponding to pure end-over-end motion (Fig. 9a), and the same object with a 10-degree “wobble” as 
described in section 4 (Fig. 9b). 




Fig. 9: Diffuse (red) and specular (yellow) intensities from a simulated RB seen from New Mexico Skies observatory. In case 
(a) no wobble (b) satellite has 10-degree wobble as defined in text. 


The most immediately evident difference between cases with and without wobble is that specular flashes occur at 
differing times. This is understandable, since specular flashes only occur when the phase angle bisector is parallel to 
a surface normal vector, so that tipping of the long axis due to wobble will change the occurrence of this condition. 



Fig. 9: Diffuse (red) and specular (yellow) intensities from a simulated RB seen from New Mexico Skies observatory. In case 
(a) no wobble (b) satellite has 10-degree wobble as defined in text. 


AZ(deg)=32 EL(deg)=7 RANGE(km)=2765 SPA(deg)=74 
object name: SL-8 R/B 04784U yr2012 day144 hr2 min39 t(sec)=900 
Group III: no wobble. dur=15min, 201 frames: based on previous SN4784 case: 6/20/12 
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If desired, the program simultaneously 
outputs a video of the RB as seen from the 
observation site. Figure 10 is an example 
of a video frame from the program’s 
output. The yellow, green, and magenta 
lines represent the directions toward the 
sun, the observer, and the angular velocity 
vector, respectively. So that the green line 
is not reduced to a point, the view is 
rendered from a location a few degrees 
away from that of the observer’s location. 

9. Discussion and conclusions 


The characteristics of realistic rotation 
states for cylindrical RBs undergoing free 
precession and forced precession have been 
examined in some detail. It is found that in 
free precession, values of the roll rate 
psidot are comparable to the rate of 
precession (phidot) of the long axis about 

the angular momentum direction in inertial space, until the object is within only a few percent of the minimum 
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Fig. 10: Output frame from Matlab observation simulation program. See text. 






energy rotation state. Furthermore, estimates of rates of internal dissipation as functions of the coning angle theta, 
or alternatively, of the energy, suggest that dissipation may slow dramatically as the minimum energy state is 
approached, so that theta may linger within several degrees of a flat spin for an extended period of time. Although 
phidot is the primary modulator of the light curve signal as flat spin is approached, it is not identical to the rotation 
rate because of the persistence of significant roll rate psidot. RBs with assymmetric mass distributions, for which 
the principal axes of inertia are offset from the symmetry axes, will be indistinguishable from RBs that are 
undergoing free precession, unless they have noticeably non-axisymmetric reflection characteristics. For Soviet SL- 
8s, the saddle tanks may produce roll-rate signatures of -25% variation in light curve intensity, which may be 
observable and could provide confirmation of estimated values for phidot and theta, through the relationship given 
in equation (3). Alternatively, if detectable in telescopic observations, colored insignias, present typically on one 
side of RBs at launch, would generate signals providing an ideal measure of phidot. The only torques large enough 
to generate significant forced precession are the gravity gradient torque and the magnetically induced eddy current 
torque. However, for typical rotation rates near flat-spin conditions, these torques can only cause precession periods 
of weeks to months, whereas free precession rates are comparable to the rotation rate. Published claims of rotation 
periods, precession periods, and coning angles that do not conform to these specifications (e.g., Toshi’s paper) are 
therefore necessarily erroneous. When slowing rotation rates approach the mean motion, large deviations of 
orientation from the commonly expected state where the long axis is directed toward the geocenter may be common, 
due to the proximity of spin-orbit commensurabilities. Numerical integrations of the full Euler equations of motion 
for an RB are under development. These integrations include the effects of the eddy currents, gravity gradient 
torque, and internal friction. 
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